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Abstract. In this paper we study sequences of matrix polynomials that sa- 
tisfy a non-symmetric recurrence relation. To study this kind of sequences we 
use a vector interpretation of the matrix orthogonality. In the context of these 
sequences of matrix polynomials we introduce the concept of the generalized 
matrix Nevai class and we give the ratio asymptotics between two consecutive 
polynomials belonging to this class. We study the generalized matrix Cheby- 
shev polynomials and we deduce its explicit expression as well as we show some 
illustrative examples. The concept of a Dirac delta functional is introduced. 
We show how the vector model that includes a Dirac delta functional is a 
representation of a discrete Sobolev inner product. It also allows to reinter- 
pret such perturbations in the usual matrix Nevai class. Finally, the relative 
asymptotics between a polynomial in the generalized matrix Nevai class and 
a polynomial that is orthogonal to a modification of the corresponding matrix 
measure by the addition of a Dirac delta functional is deduced. 



1. Introduction 

In the last decade the asymptotic behavior of matrix orthonormal polynomials, 
the distribution of their zeros as well as their connection with matrix quadrature 
formulas have paid a increasing attention by many researchers. A big effort was 
done in this direction in the framework of the analytic theory of such polynomials 
by A. J. Duran, W. Van Assche and coworkers, among others (cf. [9j H3 l H4 l [T6 | 125 ] ) . 

In this work we study outer ratio asymptotics for matrix orthogonal polynomials 
belonging to a new class, the so called generalized matrix Nevai class, and for these 
matrix orthogonal polynomials we obtain some new analytic results. 

Definition 1. Let A, B, and C be matrices, with A and C non-singular matrices 
of dimension N x N . A sequence of matrix polynomials {V rn } m £fi defined by 

zV m (z) = A m V m+1 (z) + B m V m (z) + C m V m _ 1 (z), m > 0, (1) 

with A m a non-singular lower triangular matrix and C' m a non-singular upper- 
triangular matrix, belongs to the generalized matrix Nevai class M(A, B, C) if 

lim A m = A, lim B m — B, and lim C. m = C. 

m— foo m— »oo m— foo 
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We say that a matrix of measures M belongs to the generalized matrix Nevai class 
M(A,B,C) if some of the corresponding sequences of matrix orthogonal polynomials 
belongs to M(A,B,C). 

In order to to study the generalized matrix Nevai class we will recover a vector 
interpretation of the matrix orthogonality that was presented for the first time 
in [7j. Let (P^)* be the linear space of vector linear functionals defined on the linear 
space F N of vector polynomials with complex coefficients. (P^)* is said to be the 



dual space. A vector of functionals 11 
defined by 



i N ] T acting in P N over M N xn(C) is 



U(T) := {U.y 1 ) 



(u N ,Pi) 



(u N ,p N ) 



where "." means the symbolic product of the vectors XL and T 7 ' , where 7 T — 
[pi ■ ■ ■ Pn]i Pi G the standard linear space of polynomials with complex coeffi- 
cients. The degree of 7 is given by deg(T) = |_( m & x {deg pj})/N\, where |_-J 

j=l,...,N 

represents the integer part of a real number. 
Given a polynomial h, with deg h — N, the set 

{1, x, . . . , I** -1 , h(x), xh(x), . . . , x N ~ 1 h(x), h 2 (x), xh 2 (x), . . .} 

is a basis for the linear space of polynomials, P. Furthermore, {CPjjjgN, with 7j 



defined by ^-(a;) = {h{x)yy {x), where ? {x) = [Ix 
the linear space of vector polynomials P N . 

The vector of linear functionals (x k VL) acting in 



„n-ut 



is also a basis for 



N 



{x k U)(9) := ((x k U).y T f 



over Mnxn(C) is defined by 

u(x k y) . 



So, with this definition and taking into account that is a basis for the 

linear space of vector polynomials P^, the j-th moment associated with the vector 
of linear functionals x k U is given by (x fc U)(J' J ) = U k . 
The Hankel matrices associated with U are the matrices 



D„ 



U 



1J 



11 



u 



2 m 



m G N, 



where Uj is the j-th moment associated with the vector of linear functionals U. U 
is said to be quasi- definite if all the leading principal submatrices of D m , m G N, 
are non-singular. 

A vector sequence of polynomials {23 m }„ l6 N, with degree of 23 m equal to m, is 
said to be left- orthogonal with respect to the vector of linear functionals 11 if 

(h k U) (S m ) = A„A, m , k = 0, . . . ,m - 1, m G N, 

where 5k t m is the Kronecker delta and A m is a non-singular upper triangular matrix. 
Notice that, we can always write H m in the matrix form 



2m(ar) = V m (h(x))-J> {x), 
where V m is &m degree TV x N matrix polynomial and T'q(x) = [1 x 



„N-1]T 
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Similarly, a sequence of matrix polynomials {G m }meN, with degree of G m equal to 
m, is said to be right-orthogonal with respect to the vector of linear functionals U if 

(G^{h{x))U) = e m Sj, m , j = 0, 1, . . . ,m - 1, m G N 

where O m is a non-singular lower triangular matrix. 

In [7] necessary and sufficient conditions for the quasi-definiteness of U, i.e., for 
the existence of a vector (matrix) sequence of polynomials left-orthogonal (right- 
orthogonal) with respect to the vector of linear functionals U are obtained. 

These sequences of polynomials satisfy non-symmetric three-term recurrence re- 
lations. So, if { < B m }meN i s a vector sequence of polynomials left-orthogonal with 
respect to U and if {G m } m£ N is a sequence of matrix polynomials, then there ex- 
ist sequences of numerical matrices {A m } m( zf$, {£? m } m gN, and {G m } mg N, with A m 
a non-singular lower triangular matrix and C m a non-singular upper triangular 
matrix, such that 

h(x)V m (x) = A m V m+1 (x) + B m -B m (x) + Cm'Brn-lix), TO > 1, (2) 
with = Oixat and T> Q (x) = 9 (x) , where 3> (x) = [lx ■■■ /^] T and 

zG n (z) = G„_i(z)A n _! + G n (z)B n + G n+1 (z)C n+1 , n > 1 , (3) 

with G-i{z) = Ojvxat and G (z) = U(T ) _1 - 

Notice that these three-term recurrence relations completely characterize each 
type of orthogonality. 

Furthermore, right and left vector orthogonality are connected with right and left 
matrix orthogonality. Indeed, consider the generalized Markov matrix function 2f 
associated with It defined by 



T(z) := U, 



I — h{x) 



z-h(x)) "' ( U x> z-h(x)l 



V U x) z -h{x)l \ U x ' z-h(x)l 



with z such that < \z\ for every x € L where L = Uj—i^,,^ supp vP x . Here 

1L X represents the action of U on the variable x and Tq(x) the same as before In 
fact, see [7], the matrix sequence {G ra } nS N and the vector sequence {23 m }meN are 
bi-orthogonal with respect to U, i.e., 

{{G n {h(x))) T U x ){T> m ) = I NxN 5 n , m , n,meN, 
if and only if matrix the sequences {G„}„ e N and {V m }mgM, where 

3 m W = V m (h(z))9 (z), (4) 
are bi-orthogonal with respect to if, i.e. 
1 



2ni 



V m (z)3 r (z)G n (z)dz = I NxN 5 n m , n, m e N, 



c 



where C is a closed path in {z e C : \z\ > \h(x)\,x G L}. 

The sequences of matrix polynomials {V m } m £f$ and {G m } m6 N presented here are 
orthogonal with respect to a matrix of measures which is not necessarily positive 
definite. These sequences satisfy the three-term recurrence relations (JTJ) and 
respectively. On the other hand, these recurrence relations yield a characterization 
of right and left matrix orthogonality. 
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The sequences of matrix polynomials {!Bil'} me s and {Gm } m gM given by 

»£>(*) = Ug ( W*)-WW*)) 

14 



G m+1 (z)-G m+1 (M*)) 



are said to be the sequences of associated polynomials of the first kind for U and 
{'BmjmeN and {G m } m gN, respectively. Here tC^ represents the action of 11 on the 
variable x. 

As a consequence of the definition of associated polynomials of the first kind for 
23 m and G m , we get 

'23 TO+ i(a;) > 



V m+1 (z)?(z)-%$(z)=U x ( 



z — h(x) 



(5) 
(6) 



J(z)G m+1 (z) - GW(«) = (G T m+l {h(x))U x ) (^|y 

The sequences of associated polynomials of the first kind {S™ } ra£ n, {Gm } me n 
satisfy respectively the three-term recurrence relations ([1]) and (J3]) with initial con- 
ditions S^z) = Oatxjv, $£\z) = A^ 1 and G^z) = NxN , Gg\z) = Cf 1 . 

The sequences of matrix polynomials {V m } m< =n and {G m } m gN satisfy a Chris- 
toffel-Darboux type formula 

m 

(x- z)J~] G k (z)Vk(x) = G m (z)A m V m+ i(x) - G m+ i(z)C m+ iV m (x), (7) 

fe=0 

with x, z € C and its confluent form 

m 

^Gfe^)^^) = G m (!)yl m l/: +1 (i) - G m+ i(aOC m +iVM, (8) 

k=Q 

with ieC. The Christoffel-Darboux formula characterizes the matrix orthogona- 
lity (cf. [B]) and allows us to deduce the following result. 

Theorem 1 (Liouville-Ostrogradski type formula). Let {V m } men and {G m }meN be 
the sequences of matrix polynomials bi-orthogonal with respect to £F with {V m }meN 
defined by J3J). Let {23m }me« and {Gm }meN be, respectively, the sequences of 
associated matrix polynomials of the first kind for {23 m }meN and {G m } m eN- Then, 

^G m -V m+1 G^ ) _ 1 =A- i (9) 
where A m is the non-singular coefficient that appears in the recurrence relation . 

Proof. The sequences of polynomials {V m } me ®, {G m } meN , {23^ ) _i}meN, and {G^_i}m£N 
satisfy, respectively, the recurrence relations ([1} and ([3]) with standard initial con- 
ditions. 

To prove this result we proceed by induction. For m = the result follows from 
the initial conditions. We assume the formula 

-B^G p -V p+l G { ;l x =A-\ 

is true for p = X, . . . ,m — 1. To prove that this relation is also valid for p = m we 
consider the following steps. First, we use the recurrence relation in 2$m and in 
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V m +i, i.e., 

^m^G m — V m j t iGU_ 1 — A^ (zInxn ~ B m )(H^_ 1 G rn — VmG^^-^) 

~ An ^"»(B m -2^m — Kn— l^m— l)' 

Second, we prove that CB^_iG m - V m G < ^ > _ 1 = Onxn- 

Multiplying by G m in the right side of ([5]) and multiplying by V m in the left side 
of (J6j> , changing m by m — 1 in the relations (0 and (JOJ) , and then subtracting these 
equations, we get 

Then, 

%%_ 1 ( z )G m (z)-V m (z)GW_ 1 (z) 

= V m (z)(Gl(h(x))U x ) (M^-U* (^gy) G m (z). 

Adding and subtracting (Gj n (h(x))U x ) ( -^zj^j ) m the last relation and taking in 
consideration the left and right orthogonalities, the result follows. 
Using the above result we have 

"B^Gm — V m +lG^_ l = —A^G'm (^B^ m L 2 G m — Vm-lG^L^J . (10) 

Again, using the recurrence relations for G m and G^ > _ 1 in r B^_ 2 G m —V m -iG^_ 1 
we get 

2 r ^_ 2 G f m - Vm-lGm-1 = (^-2^-1 - V m -iG^_ 2 ){zI N xN - B m -l)C^ 

+ (V m -iG^_ 3 - H < ^ ) _ 2 G m -2)A m -2C m 1 . 
Since S^_ 2 G m _i - V m -iG£_ 2 = NxN , we get 

Using this relation in (ITU)) we obtain 

25^G m - Um+iG^Li = -A m 1 G m (U m _iG| r 1 l L 3 - ®^_ 2 G m _2)A m _2G TO 1 . 
According to the hypothesis of induction the result follows. □ 

Let U be a quasi-definite vector of linear functionals and let {V m } m eN and 
{G m } me N be sequences of bi-orthogonal polynomials. We denote the kernel poly- 
nomial by 

m— 1 

K m (x,y) = Y,Gk(h(3,))V k {h{x)). 

k=0 

Notice that, even though we don't have K m (x,y) — K^y^x) like in the matrix 
symmetric case, the reproducing property for the kernel holds. 
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Theorem 2. Let U be a quasi- definite vector of linear functionals, {¥> m } 

mdN ana 

{G m }meN be, respectively, the left vector and right matrix orthogonal polynomials 
with respect to 11. Given a vector polynomial ir £ P N of degree m, i.e., 



(11) 



fe=0 



thenir(x) = (K% +1 (x, z)U z )(tt(z))%{x). 

Proof. From (fTT|) , using the bi-orthogonality and taking into account L n = G^ll, 
we get ft? = L k (ir(z)) = (C%(h(z))U x )(ir(z)). Then, 



n(x) =J2(GZ(h(z))U z )(ir(z))'B k (x) 



k=0 



But, 



(Gl(h(z))U z )(Tr(z))-B k (x) = ((G k (h(z))V k (h(x))) T U z ))(7r(z))7 (x). 



Hence, the result follows. 

Let J be the block matrix defined by 



□ 





'B 


A 







Ci 


Bi 


A\ 


,/ = 













c 2 


B 2 



(12) 



that is known in the literature as N-Jacobi matrix. When the matrix polynomials 
satisfy a symmetric recurrence relation, it was proved in [15] that the zeros of the 
m-th orthogonal polynomial are the eigenvalues of J m . This result can be general- 
ized for sequences of orthogonal polynomials satisfying non-symmetric recurrence 
relations. Thus, for m G N, the zeros of the matrix polynomials G rn and V m arc 
the zeros of the polynomial det(tI m N Xm N — J m ), with the same multiplicity, where 
IrnNxmN is the identity matrix of dimension mN x mN and J m is the truncated 
TV-block Jacobi matrix of dimension mN x mN. 

Also, for a matrix of measures W and for any polynomial V with degree less 
than or equal to 2m — 1 the following quadrature formula holds 



f V(h(x))dW(h(x)) =J2V{x m , k )T m>k 
J fe=i 



where zeros of the matrix polynomial V m , in general, 

complex numbers as well as s < mN, and T m ^ k are the matrices 

h 



■ , i | , t - ,„)( ^Adj{V m (x)))^- l >{x m , k )^i(^k), 

(det (VwWjj^^im,*) 

for fc = 1, . . . , s, where l k is the multiplicity of the zero x m . k . 

Using the above quadrature formula, in [7] we have obtained the following asym- 
ptotic result: 

lim V-^z)^ t {z) = Hz) 



locally uniformly in C \ T, where T — Dn^qMn, Mn = U„>j\r{zeros of V n } . 
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The structure of the manuscript is as follows. In section [21 sequences of matrix 
orthogonal polynomials belonging to the generalized matrix Nevai class are studied. 
Furthermore, the outer ratio asymptotics of two consecutive polynomials belonging 
to this class is obtained. We also study the generalized matrix Chebyshev polynomi- 
als and we present their explicit formulas as well as the corresponding generalized 
Markov function. The example presented in this section does not belong to the 
cases studied by A. J. Duran in [12] and also can't be converted in these cases by 
the method presented by H. Dette and coworkers (see [5]). 

In section [31 we introduce a modification of a vector of linear functionals by 
adding a Dirac delta. This yields a reinterpretation and an extension of a pertur- 
bation in the usual matrix Nevai class (see f24j [27l [28]). On the other hand, the 
meaning of any discrete Sobolev inner product in vector terms is clarified and this 
is clearly related with sequences of polynomials satisfying higher order recurrence 
relations (cf. [3 IB H31 [231 [301 13U 132 ) ■ We also find necessary and sufficient con- 
ditions for the quasi-definiteness of the modified functional. The quasi-definiteness 
conditions obtained by the authors in p] [2l [21] for some special examples coin- 
cide with our general result. To conclude this section, we describe the generalized 
Markov function associated to a modification by a Dirac delta functional. 

Finally, in section 2] we present the relative asymptotics between a sequence 
of matrix polynomials that belongs to the generalized matrix Nevai class and a 
sequence of matrix polynomials that is orthogonal with respect to a modification 
by a Dirac delta functional of this class. This result generalizes those obtained by 
F. Marcellan and coworkers (see [27]). 



2. Generalized matrix Nevai class 

Let M be a matrix of measures in the generalized matrix Nevai class M(A, B, C). 
Notice that M can belong to several Nevai classes because of the non-uniqueness of 
the corresponding sequences of orthogonal polynomials. 

If A and C are non-singular matrices we can introduce the sequence of matrix 
polynomials {U^ B ' C } m< z^ defined by the recurrence formula 

zU£ B - c (z) = Auiff{z) + BU^ c (z) + CU^fiz), m > 1, (13) 

with initial conditions U ' B,c (z) — Inxn and U_{ B ' C (z) = 0n x n- According to an 
extension of the matrix Favard's theorem (see [7]) this sequence is orthogonal with 
respect to a matrix of measures ~Ma.b,c that is not necessarily positive definite. 
This sequence of matrix polynomials is said to be the sequence of matrix generalized 
second kind Chebyshev polynomials. 

The continued fraction associated with (|13j) or, equivalently, the generalized 
Markov function 3 a, b,c is given by 

3a,b,c{ z ) = J— r m * 1 n 

(ZlNxN — a) — A j O 

zI NxN -B- A— C 

ZlNxN — 

where 1/X denotes the inverse of the matrix X. 

Matrix continued fractions of different types were studied by many authors (cf. [31 
[T51 121)1 [33]). In [33] the reader can find a detailed study about matrix continued 
fractions and matrix Chebyshev polynomials, where the author emphasyzes how 
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continued fractions are used to develop the notion of matrix Chebyshev polynomials 
in some symmetric cases. 

Now, let us consider a sequence of vector polynomials {23 m }meN left-orthogonal 
with respect to the vector of linear functionals U satisfying the recurrence relation 

h{z)T m {z) = A% n+1 {z) + B-B m {z) + CV m -i(z), m > 1, 

with initial conditions CB_i(z) = Oat x i and T>o{z) = ^o(z), where A and C are 
non-singular matrices. It is straightforward to prove that the sequence of ma- 
trix polynomials {U m }meN, defined by H m (z) = V m {h{z))'3 > o{z), and the sequence 
{23m'}meN, of associated polynomials of the first kind for II and {23 m } me N, satisfies 
the same recurrence relation with the following initial conditions V-i(z) = Onxn, 
Vo(z) — Inxn, 'S'-l(z) = Onxn, and 1}q~\z) = A^ 1 , respectively. These conditions 
are said to be the standard ones. 

Rewriting the recurrence equations for {KjmeN and {B„ } me g in a blocks 
matrix form, we have 



AV m+1 



AT, 



(i)' 



3 



(i) 



zInxn — B 
Inxn 



-C 

OjVxIV 



V m 

V m -1 



(1) 
m — 

(i) 



Since the matrix A is non-singular the last equation is equivalent to 



V m +i 

V m 



(i) 



3 



(i) 



m— 1 

Writing the last equation as L 



A x (zI N k N 
Inxn 



B) 



-A~ X C 
Onxn 



V„. 
V m - 



(i) 

m — 

(i) 



TL m _i where 



Lm — 



Vm+1 



and T = 



A- l (zI Ny . N -B) 
Inxn 



-A~ X C 

OnxN 



we have L m = T m L with 



A~ 



1 (zIn x n 
Inxn 



B) 



A- 1 ' 

OjVxAT 



For some particular choices of A, B, and C the matrix T' m has the spectral decom- 
position T' m = SD m S~ 1 where D is a diagonal. 

Using this decomposition we can determine L m and then, obtain V m and 3„Li- 
By a straightforward calculation we obtain V^ 1 and then taking the limit 

when m — > oo, we get the generalized Markov function Ja.b.c- 

As a sake of example, in the case of matrix polynomials of dimension 2 x 2, if 
we consider the same problem but with different initial conditions, for instance, 
V {z) = P, Viz) = M + Qz, %W(z) = 2x2 , and %g\z) = Q, with 



p = 


Pll 


Pl2 


, M = 


mu 


m 12 




P21 


P22_ 




m 2 i 


m 2 2_ 



and Q = 



In 

921 



912 
922 



we can relate the generalized Markov function associated with these new initial 
conditions, 9\4,b,c, with ^a.b.c, in the following way: 
- The block matrix L m is given by 

L m = T Lq = T Lq = T m LQL 1 Lq = L m L Lq. 
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Using this relation we get, 

V m (z) = V m {z)P + 'B^ ) _ 1 (z)[AM + BP+(AQ-P)z} 



3a,b,c{z) = P- 1 y A ,Bfi{z)AQ + M- 1 Q + p- 1 B- 1 AQ + {I - P- 1 AQ)z. 

In the next example, we illustrate how can we determine the generalized Markov 
function associated with generalized second kind Chebyshev polynomials for a par- 
ticular choice of the recurrence coefficients with standard initial conditions. 



Example 1. Let us consider A, B, and C as 

-1 



A = 



1 
1 



B 



1 





-1 



and C 



The matrices T and Lq are 
"1 + 2 



2 + 1 



1 



-1 
1 





and L 



1 + 2 
-1 
1 





-1 
1 



1 

1 + 2 1 



1 



Matrix T have the following eigenvalues 
1 



Ai = -(1 + 2 - v/(* + l) 2 -4), A 2 = -(1 + 2 + y/( z + i)2-4), 

A 3 = i(l + 2- v/(2 + l) 2 + 4), A 4 = i(l + 2 +v /(2 + l) 2 +4) 
and the corresponding eigenvectors are 



vi 



V2 



1 



2, -(1 + 2- v/(2 + l) 2 +4), 



-1-2+^/(2 + 1)2+4'" 



-2,i(l + 2 + v/(2 + l) 2 +4),- 



1 + 2+^/(2 + 1)2+4 



• 1 



V3 



0,-(l + z- y/(z + l) 2 -4), 0,1 
0,^(1 + 2+ ^(z + l) 2 -4), 0,1 



Then, the matrix T' m has the following spectral decomposition 

T m = SD m S~ 1 , 

where D = diagonal [Ai, A 2 , A 3 , A 4 ], for m e N, and S = \v-\\ v%\ v$\ v^], where \ 
and Vi, for i = 1, . . . , 4 are, respectively, the eigenvalues and the eigenvectors of the 
matrix D. Thus we can determine L m+ \ and then, obtain V m and r B^_ 1 as follows 



Vm(z) = -^E m+2 (z) 





1 + 2 



I- E m+ i(z) 
+ F m+1 (z) 



1 + 2 

-\ q 





i(2 + 2)2 1 + 2 



+ E m (z) 
-F m {z) 



1 




1 + 2 1 
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where 

E m {z) 



E, 



m+2 



+ E m+ i(z) 



1 




+ F m+1 (z) 





i(l + z) 1 



,(!+£+ + I) 2 + 4) m - (1 + z - y/(z + l) 2 + 4) T 
^(l + z) 2 + 4 



= 2" 



(1 + z + ^(z + l) 2 - 4)" 1 - (1 + z - v /(z + l) 2 -4) ? 



By a straightforward calculation we obtain V m lr B^ > _ 1 . 
m — > oo, we get 



Taking the limit when 



3a,b,c{z) = 



l+z+V(l+z) 2 +4 



4+(l+2-^(l + z) 2 +4)(l + z- A /(l + z) 2 -4- % /(l + z) 2 +4) 
(l+2 + ^/(l+2) 2 -4)(l + z+^/(l + 2) 2 +4) 



l + z + v /( 1 + z ) 2 - 4 



This example inspired us to prove the ratio asymptotics between two consecu- 
tive polynomials in the generalized matrix Nevai class. Indeed, without loss of 
generality, we consider a sequence of matrix polynomials {V m } m ^ defined by ([T]) 
such that Vq(z) = Inxn- But first, we need to introduce the following result. 

Lemma 1. Let {V m } m ^fs be a sequence of matrix orthogonal polynomials in the gen- 
eralized matrix Nevai class M(A,B,C). Then, there exists a positive constant M, 
which does not depend on m, such that their zeros x m ^ are contained in a disk 
D = {z e C : \z\ < M}. 

Proof. Let us consider the TV-block Jacobi matrix, J (see (|12p). associated with the 
recurrence relation ([T]). Remember also, that the zeros of V m are the eigenvalues 
of J rn where J m is the truncated matrix of J, with dimension mN x mN. 

Taking into account that the sequences {A m } m£ fi, {-B m } m gN, {C m } m eN con- 
verge, and using the Gershgorin disk theorem for the location of eigenvalues, it 
follows that there exists M > such that if x„ l: k is a zero of V m then x m> k € D 
where D = {z € C : \z\ < M}. So, L defined by ' 



T = n N > M N , M N = U m >jv{zcros of V m }, 
is contained in D and supp(M / ) cfcf. 



□ 



Theorem 3. Let {V m }meN be a sequence of matrix polynomials left- orthogonal 
with respect to the matrix of measures M and satisfying the three-term recurrence 
relation (fT]). Assume that lim A m — A, lim B m = B, and lim C m = C with 

A and C non-singular matrices. Then, 



lim V m _ 1 {z)V- 1 {z)A^_ l = 5 a ,b,c{z), zeC\T, 



(14) 



where 3a,b,c is the Markov transform of the matrix of measures for the generalized 
second kind Chebyshev polynomials. Moreover, the convergence is locally uniformly 
for compact subsets ofC\T, where T — C\N>oMpf, M/v = U m >NZ m , and Z m is the 
set of the zeros of V m . 



RELATIVE ASYMPTOTICS FOR ORTHOGONAL MATRIX POLYNOMIALS 



11 



Proof. First, we consider the sequence of discrete matrix measures 
{fJ-m}m£N defined by 

s 

Mm = / ] Vm-l(h(ym,k))Fm,kGm-l(h(ym,k)) &y m ,k-> m > 0, 
k=l 

where y m ,k are complex numbers such that h(y mi k) — %m,ki with x mi k, k = 1, . . . , s, 
the zeros of the polynomial V m , and the matrix T mt k is given by 

r / fc (Adj (V m (x)))^- 1 Hx m , k )'B { J l ) _ 1 (x m ^ 
m ' k ~ (det(V m (x)))^)(x m , k ) i! 

Ik being the multiplicity of the zero x m< k, h < N, and {23^Li}meN the sequence 
of associated polynomials of the first kind for {23 m } m(E N and It. Notice that T m k 
is the weight in the quadrature formula. Hence, it follows that 



d/j, m (h(x)) = InxN, for m > 
The decomposition, (cf. [!§]). 

V m ^{z)V-\z) = VC m , fc 

with 



K — U 



r = ?fcKr 1 -l(x m ,fc)(Adj(T/ m ft))) ( ' fc ~ 1) (^m,fc) 

m - k (det(V m (t)))^)(x m , k ) 
is always possible even though the zeros of V m are complex or have multiplicity 
greater than one (see [5J [7J [H]). Then, we have 

! _ ? fc y m -i(x m , fc )(Adj (V m (t)))(' fc - 1 Ha: m ,fc)^ 1 _ 1 

m ' fc m - 1_ (det^i)))^)^) 

Applying the generalized Liouville-Ostrogradski formula ([9]) and taking into account 
that for every b, a zero of a matrix polynomial V m , (see [10] ) 

V m (b) (Adj (K^)))^ 1 ' (b) = (Adj (V TO (t))) (p - 1) (b)V m (b) = WxJV , 
we obtain 

Cm,fc^TO-l = ^ra-l(^(2/m,fc))r TO ,fcG m -l(/l(2/ m ,)fc)). 

From the definition of the matrix of measures fj, m we get 

v m ^{z)v m 1 {z)A m \ = .ec\r. 



Let us consider the generalized Chebyshev matrix polynomials of second kind 
{Um' B ' C }meN defined by (TT3]). We can prove by induction that 

lim / Uf< B ' C (h(t))dn m (h(t)) = I f™' lZ \ Z n (15) 



m->oo y * v v " ^ mv v ;; [ OatxAT, for Z ^ 

To do it, we just use the same technicalities as in [TJ]. 
We are now ready to prove 

m->oo J z — h[X) 
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> C > 0, I > 0, (16) 

2 



If not, we can find a complex number z £ C \ T, an increasing sequence of nonneg- 
ative integers {rt;}; 6 N, and a positive constant C such that 

f dfi ni (h(x)) 

J z~h(x) "^' cW 

where |j . H2 denotes the spectral norm of a matrix, i.e., 

||A|| 2 = max {VA: A is a eigenvalue of A* A}. 

Since {/x m } m£ N is a sequence of matrices of measures with support contained in 
a disk D (see lemma [l} and taking into account that J dfi m = Inxn, by using 
the Banach-Alaoglu theorem, we can obtain a subsequence {r;}; 6 N from {n;}z e N, 
defined on a curve 7a/ contained in the disk D, with the same fc-th moments of the 
vector of linear functionals IX, for k < 2r; — 1, such that 

lim / f(h(x))dfi n (h(x)) = ^ [ f(h(z))U x ( )dz, 
l-xx, J lM 2ni J JM \z-h(x)J 

for any continuous matrix function / defined in D. 
Hence, taking f(h(x)) = U A > B ' c (h(x)), we have 

A,B,C lu i u , /, / n _ 1 f TT A,B,C( h ( \ \ rw ( ^o( x ) 



Km / uC'^(h(x))d^ ri (h( x )) = — U^ D ^(h(z))U 
From (IT5|l we have 

1 /" tf, A ' B - c (M^)ic/ To(x) ^ - f lNxN ' for/ = °' 



27ri7 7M « v y " x \z-h(x)J \ Onxn, for/ 7^0. 

But, the sequence of matrix polynomials {U^' B }meN is orthogonal with respect 
to ^a.b.c- Since {C/^ ,B,c } m gi is a basis of the linear space of matrix polynomials 
we get that (fl6|) is not possible. Each entry of the matrix J d ^l^^ is uniformly 
bounded on compact sets of C \ T. Then, according to the Stieltjes-Vitali theorem, 
we get the uniform convergence. □ 

Notice that we have analogous results of Lemma []] and Theorem [3] for the se- 
quence of matrix polynomials {G TO } m gN- 

Corollary 1. Under the hypothesis of Theorem\^we have that 

limjV^zW-Hz^A-l, = 5% tC {z) (17) 

on compact subsets ofC\T, for k = 1, 2, . . . . 

The locally uniformly convergence in (|17p means that every entry of the left 
hand-side of (|17p is locally uniformly convergent to its corresponding entry in the 
right hand-side of ([TT]) . 

3. Delta functionals 

In this section we deal with a vector of linear functionals that results of a modi- 
fication by a Dirac delta functional and we will illustrate it with a nice application 
related to Sobolev inner products. Thus, we can reinterpret these inner products 
in a vectorial form and it is a motivation for the study of these modifications. 

In a generic way, we can say that a Dirac delta functional is a vector of linear 
functionals where the functional components are linear combinations of Dirac deltas 
and their derivatives in a finite set of points in the real line. 



RELATIVE ASYMPTOTICS FOR ORTHOGONAL MATRIX POLYNOMIALS 



13 



If the polynomial h of fixed degree N consider in the previous sections is such that 

M 

h(x) = l[(x-c j ) M *+\ (18) 

3=1 

where Mj + 1 is the multiplicity of each Cj € N as a zero of h, then we can define 
a new vector of linear functionals as the result of a modification by the addition of 
a Dirac delta functional with respect to h as follows. 

Definition 2. The vector of linear functionals li defined by 

U = U + A5, (19) 

with A a numerical matrix of dimension N x N where 

6 -IS S' ■■■ 8 {Ml) 8 5' ■■■ 5 {M2) ■■■ 6 5' ■■■ 6 {Mm) ] t 

where N — M + Xwii * s ca ^d vector of linear functionals modified by a Dirac 
delta functional associated with h. 

First, as a motivation for the study of these modifications, we will consider the 
following example. It gives a vectorial reinterpretation of a Sobolev inner product. 
It is important to refer that with the same technics presented in this example, a 
vector reinterpretation for any general discrete Sobolev inner product holds. 

Example 2 (see [21]). Let us consider the discrete Sobolev inner product 

(f,g) s := J fgdn+\f'(pW(0), where X e R + . (20) 

To establish the parallelism between vector orthogonality and Sobolev inner products, 
let {p n }n£N be a sequence of scalar polynomials orthonormal with respect to the 
inner product (|20|) . i.e. 

(Pn,x k ) s = 0, fc = 0,...,n-l, (21) 
(Pn,x n ) s + 0,neff. 

Notice that the multiplication by the polynomial h{x) — x 2 yields a symmetric 
operator with respect to the above Sobolev inner product, i.e., 

(x 2 f,g) s = (f : x 2 g) s ,\/f,geP, 

and, as a consequence, the corresponding sequence of orthonormal polynomials 
{Pn}nen satisfies, for n > 0, a five-term recurrence relation 

X 2 Pn+l = C n+3 ,2Pn+3 + C„ +2 ,lPn+2 + Cn+lfiPn+l + C n +l,lPn + Cn+lflPn-l ■ (22) 

On the other hand, taking into account the vectorial approach given in [7j, the 
recurrence relation (|22[) yields the following vector expression 

x H m (x) = Am+i^m+i ( x ) + B m T> m (x) + Al l , B m _ 1 (x),m>0, 

where ~B m (x) = [p 2m (x) P2m+i(x)] T , 



A m — 



C2m,2 
C2m,l C2m+1,2 



id B m = 



C2m,0 C2m+l,l 
C2m+l,l C2m+l,o 



Notice that B m — Bj n as well as we are dealing with a symmetric case. Then, 
by Favard's type theorem (see ^\), there exists a vector of linear functionals 11 = 
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[u 1 u 2 ] T such that the sequence of vector polynomials {25 m } m gN is left- orthogonal 
with respect to IX and 



(x 2k U) (s m ) - 2x2 , fc = 0,l,...,m-l, 
(x 2m U) (% m ) = A m , m e N, 



(23) 



where A m is a non-singular upper triangular matrix. 

From the definition of a vector of linear functionals, the previous orthogonality 
conditions are given explicitly, for all k — 0, ... , to — 1 and to €E N by 



(x 2k U) (E 
(x 2m U) (S m ) 



(u 1 ,x 2k p 2m ) 

(u 1 , X 2k p 2m +l) 



(u\x 2m p 2m ) 
(u\x 2m p 2m+1 ) 



(u 2 ,X 2k p 2m ) 
(u 2 , X 2k p2,n+l) _ 

i^,X 2m p 2 m) 

(u 2 ,x 2m p 2m +i) 








0" 












• 


• 







• 



The vector of linear functionals U represents a Sobolev inner product like J20|) 
only if it is a modification by a Dirac delta functional, as we will describe in the 
sequel. In fact, the vector of linear functionals U has the representation 



U = U + AS 



(24) 



where 



U 



A 





-A 



[So S' ] T 



and where u is a linear functional on the linear space of scalar polynomials P 
such that u(p(x)) = fjp(x)d[J,(x). Here fi is the weight function that appears in 
the Sobolev discrete inner product (|20|) . 

To illustrate that the orthogonality conditions (|2ip and (|23p are equivalent, first, 
take k = in (l2"3l). 



Furthermore, for all k = 1 
(x 2k U) (B 



(l,P2m)s 


{x,p 2m )s 




"0 


0" 


_(l,P2m+l)s 


{X,p 2m +l)s_ 











, m, m G N 



{x 2k ,p 2m ) s 
x 2k ,p 2m +i)s 



{x 2k+ \p2m)s 
(x 2k+1 ,p 2m+ i)s 









7^ 



Thus, these relations show us that the orthogonality conditions (I2ip and (|23p are 
equivalent and it means that the discrete Sobolev inner product (|20[) can &e repre- 
sented in the vectorial form by 



Motivated by this example the first question naturally imposed consists to know 
when IX defined by (|19|) is quasi-definite, i.e., when there exists a sequence of vector 
polynomials {23 m } m gN left-orthogonal with respect to IX. Then, the next step is to 
obtain necessary and sufficient conditions for the quasi-definiteness of IX. 

In the sequel, we denote by {23 m } mS N the vector sequence of polynomials left- 
orthogonal with respect to IX and {23 m } mS N the sequence of polynomials associated 
with IX defined by (TIT))) , i.e., the sequence of polynomials satisfying 



IX S m = 



JNxN 



, TO > 1, 



and IX ^23o^ is a non-singular matrix. 



(25) 
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As above, we will denote by {V m }meN the sequence of matrix polynomials de- 
fined by 

B w (:c) = V m (h(x))9 (x), 
with %(x) = [la;... x N - 1 ] T . 

Before the statement of conditions that give us quasi-definiteness for U, we need 
the following auxiliary result. 

Lemma 2. Let {23 m } me N be a vector sequence of polynomials left- orthogonal with 
respect to the vector of linear functionals li and {25 m } mg N be a vector sequence 
of polynomials associated with li defined by (|19p and verifying (|25p . Then, the 
following statements hold 

(h k U)CB m ) = (h k U)(-B m ) = NxN ,k>l, 
(h k U)(B m ) = (/i fc lt)(S m ),fc>l,meN. 

Proof. The proof of this result is straightforward taking into account that for any 
vector polynomial 5 € P^, we have 

(h k S)(9) = Onxn, for k > 1, 

with 

6 = [S C1 ^ ■ • • S C2 5> C2 ... ■ ■ ■ S CM 6> CM ■ ■ ■ S^f , 
and the result follows. □ 

In the literature about discrete Sobolev orthogonal polynomials (see for exam- 
ple [J [2j gj OH [20J HU [22]) the reader could find the same conditions that we will 
achieve in the next result. Another interesting work related with this topic is [27] . 
There, the authors studied the outer relative asymptotics between a matrix polyno- 
mial belonging to a matrix Nevai class and a matrix polynomial that is orthogonal 
to a perturbation of such a matrix of measures in the matrix Nevai class. There 
the existence of matrix orthogonal polynomials with respect to perturbations in 
the matrix Nevai class is assumed. The next theorem gives necessary and sufficient 
conditions for the existence of such matrix polynomials in a more general case. 

Theorem 4. Let li be a quasi- definite vector of linear functionals, {23 m }meN be 
the vector sequence of polynomials left- orthogonal with respect to U and {23 m } m6 N 
be the associated vector sequence of polynomials, verifying (|25p . to the vector linear 
functional U defined by (|19p . The vector of linear functionals li is quasi- definite if 
and only if 

^(y (z))A T ^ m+1 ( Q , C/ ), m> 0, 
is a non-singular matrix for all c\ zero of h defined by (|18p . 

Proof. For the vector sequences {23 m } m eN and {S m } rae n, we can consider the 
sequences of matrix polynomials {V m } m eN and {V m } me iq such that 

■B m {z) = V m (h(z)p Q (z) and % m (z) = V m {h{z))9 (z), 

with CPo(^) = [1 z ■ ■ ■ z N ~ 1 ] T . Using the reproducing property for the kernel we get 

V m (h(x)) = (Kl +1 (x,z)U z )(V m (z)). 
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Similarly, we can write 

m 

V m (h(xj) = (Kl +1 ( x> z)U z )(3 m (z)) -'£{(%{h{z))A6 z )(p rn {z))V j {h{x)), 

where <5 Z is the Dirac delta functional acting on z. Notice that in the right hand 
side of the above identity 

(Gj(h(z))AS z )(V m (h(z))9 (z)) = {tif(h(ci))AS z )(y m (h(z))9o(z)), 

holds for any zero Ci of the polynomial h, independently of its multiplicity. From 
now on, we write ft(cj) = 0. 
Again, since 

8 z (V m (h(z))9 (z)) = V m (0)S z (T o (z)), 

then 

m m 

^2(Gj(h(z))AS z )^ m (z))Vj(h(x)) = K m (0) ^^(ToW)A T G,(0)^-(/i(x)). 

3=0 j=0 

Now, to analyze (Km +1 (x, z)VL z )('B m (z)) we use the definition of the kernel 
K m (x, z) = EJlo 1 Gj{h(z))Vj(h(x)) and consider the following representation for 
G m , G m {h{x)) = YJLvPTiK*)?, Pk G 3VtjVxw(C). Then, we have 

m j 

(Kl +1 (x,z)U z )CB m (z)) = ^[(^(/ l ( a ;)) fc lX 2 )(a3 m (z))^]VS(/ l (x)). 
Hence, 

(JCfi(i,4)(iW) = A m /3™F m (^)) 

if and only if {!B m }meH is left-orthogonal with respect to U z , satisfying the ortho- 
gonality conditions 

{h k U)(v„)j = A m 6 k , m , fc = 0,l,...,m, meN, 

where A m is a non-singular matrix. Thus, 

V m (h(x)) = D m V m (h(x)) - V m {Q)8 z {V {z))A 7 'K m+1 {x,ci) (26) 

where D m = A m /3™ is a non-singular matrix. Taking a; = c\ we get 

^m(0)(/ivxjv +5 z {7 n (z))K T K m+1 {c hCl )) = D m V m (0). 

{Vm}meN, is completely determined by the data if and only if 

Inxn +8 z {y {z))A T K m+1 (ci,ci),Tn > 0, 

is a non-singular matrix for any zero q of the polynomial h, independently of its 
multiplicity. □ 

If we apply this result to the previous example, then it is straightforward to 
deduce that these conditions are the same as those obtained in [3T] for the exis- 
tence of a sequence of orthogonal polynomials with respect to the Sobolev inner 
product 

In the sequel h is a polynomial of fixed degree N defined by ([TBI and sequences 
of matrix polynomials {V m }meN and {G m } meN are given by V m (z) = YlT=o a T zj 
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and Gm{z) — J^JLo^T 2 ''' These sequences are bi-orthogonal with respect to the 
generalized Markov function J and satisfy the recurrence relations ([T]) and ([3]), 
respectively. 

Theorem 5. Let 11 be a quasi- definite vector functional, 11 quasi-definite vector 
functional defined by (j 19[) . and h be a polynomial of fixed degree N. Let 1 and J 
be the generalized Markov functions associated with IX and 11, respectively. Then, 

zJ{z) = z3{z) + S x (y a {x))A T , (27) 
with y (x) = [1 x ■ ■ ■ x 1 *- 1 } 7 * and 

S=\S 5' ■■■ 5 {Ml} 5 6' ■■■ 5 (M2) ■■■ 5 6' ■■■ S ( - Mm >] t 
with N = M + J2f =1 Mj . 
Proof. From the definition of J we have 

U x {h{x)) k V Q {x) 



zJ{z) = zJ2- 



yk+l 



= z H 



k=0 

U x ((h(x)) k y (x)) + (AS x )((h(x)) k T (x)) 



k=0 

= zJ(z) + {A6 x )(7 (x)) 
as we wanted to show. □ 

On the other hand, we can also consider the sequences {V m }meN ancl {G m } m gN 
where V m {z) — J^j^o^T^ anc ^ G m (z) = Y^Lo^T^ are bi-orthogonal sequences 
with respect to the functional 5F defined by ([27| . 

The next result states the relation between {V m } m ^ and {V m } m ^N- 

Theorem 6. Let {V m } m <£N be a sequence of matrix polynomials left- orthogonal 
with respect to the generalized Markov function J and satisfying the three-term 
recurrence relation 

h(x)V m {z) = A m V m +l(z) + B m V m (z) + CmVm-l(z), TO > 1. 

Let {V m } me fi be the sequence of matrix polynomials left- orthogonal with respect to 1 
defined by (|27[) . Then, the sequence of matrix polynomials 

{V m }meN is defined by 

h(x)V m {z) = a x m+1 V m+ \(z) + a 2 m V m {z) + a^^Vm-iiz) . (28) 

with al n+1 , a^, and a^ n _ 1 such that 

^m+l = DmAjn L m ^Q A m 
a m — D m B m + L m [3™^ C m _|_i 

^m — 1 ^mCmi 



where L m = D m V m (0)[L NxN +S z {^{z))K T K m+1 {c h ci)}- 1 S z {Vq{z))A t and D r , 
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Proof. From the proof of Theorem [4] we have 

V m {h(x)) =D m V m (h{x)) 

- A»14 l (0)[/ NX Ar+^(yo(^))A T if m+ i(c ; , C/ )]- 1 ^(y (z))A T ^ m+1 ( a ;,Q). 

For the sake of simplicity, we denote by L m the matrix 

L m = D m V m (0)[I NxN +S z (Mz))A T K m+1 (c l ,c l )]- 1 S z {y {z))A T . 

Then 

V m (h(x)) = D m V m {h{x)) - L m K m+1 (x,ci). 

Using the definition of the kernel as well as the three-term recurrence relation that 
{Vm}m&i satisfies, then 

zV m {z) = [D m A m - L m G m (0)A m ]V m +i(z) 

+ [D m B m + L m G m +i(0)C m +i]V m (z) + D m C m V m -i(z). 

Then the comparison of the coefficients in ((28)) leads to the representation of a m+ i, 
a? m , and a m _ v □ 

Example 3. In the last theorem if we consider the generalized second kind Cheby- 
shev polynomials, {U^' B ' } m eN> left- orthogonal with respect to the generalized Markov 
function Ja.b.c an d satisfying the three-term recurrence relation 

zUt Bfi {z) = AU%*f{z) + BU£ B > c (z) + CU^ B f{z), m > 1, 

and consider the matrix polynomials {U^' B ' C } m ^ left-orthogonal with respect to 
3a : b,c given by 

z?a,bM z ) = z^ B ^ c (z)+S x CP (x))A T 
with y (x) = [lx ■■■ x 1 ^' 1 } 7 and 

S-\5 5' ■■■ <S (A/l) S 5' ■■■ 8 {M2) ■■■ 5 5' ■■■ S ( - Mm) ] t 
we get that 

7tA,B,Ci \ 1 TJ A,B,Cr \ , 2 TT A,B,C/ \ , 3 TT A,B,C i \ 

zU m (z) = a m+l U m ' + { (z)+a m U m > (z) + a m _ 1 J7 m L 1 (z), 
with a m+l! a m) and af n _ 1 given by 

a ln+l = 0Z)~ 1 f3™ A - LmPo'A, 
a m-l — (Pm) 1 fimC 1 

and 

L m = 0^r 1 ^ l U^ B - c (O)[I NxN +6 z (y (z))A T K m+1 (c l ,c l )}- 1 6 z (y (z))A T . 



RELATIVE ASYMPTOTICS FOR ORTHOGONAL MATRIX POLYNOMIALS 



19 



4. Relative Asymptotics 

The works about perturbations in the matrix Nevai class (see [23 EH] ) motivated 
us to study the outer relative asymptotics in the case of the generalized matrix Nevai 
class. 

First, we present how can we reinterpret the model studied by the authors from 
a vectorial point of view. Second, the outer relative asymptotics is deduced. 
Let {G„}„ g N and {23 m } me N be bi-orthogonal with respect to U, i.e., 

Gl{h{ X ))U(% m ) = I N XNSn,m, (29) 

where IX = IX + AS with A a numerical matrix of dimension N x AT, 
S=\S 5' ■■■8 (Ml) 5 5' ■ ■ ■ 8 (M2) ■ ■ ■ S 5' ■ ■ ■ S^ Mm} ] t 

° ~~ L°Cl °C1 °C1 °C 2 U c 2 U C 2 U CM U C M U C M J ' 



and h is defined by (|T8|) . 

By a straightforward calculation in (l29|) we have 



(Gl(h(x))(U + A* a ))(S TO ) = (Gl(h(x))U)(B m ) + ((G n (h( Cl ))) T Ad x )^ m ), 

where q is any zero of the polynomial h. 
Since h(c{) — and 

((G„(0)) T A^)(S m ) = V m (0)S x (y o )A T G n (Q). 

Thus 

(Gl(h(x))(U + Atf x ))(S m ) = (G£(/i(z))U)(5 m ) + K n (0)^(T o )A T G n (0). 
By the bi-orthogonality, we have the following matrix interpretation 



1 

2^1 



V m {z)1(z)G m {z)dz 



= [ V m (z)3(z)G m (z)dz + V m (0)6 z (J> (z))A T G n {0) 

where we get the same model studied by F. Marcellan, M. Pihar, and H. O. Yakhlef 
when we choose all zeros of the polynomial h equal to c, i.e., h(x) = [x — c) . 

Next, we will deduce the outer relative asymptotics of {V m (V m )^ 1 } me ^ when 
{VmlmgM belongs to the generalized matrix Nevai class. To prove it, the following 
lemma is needed as an auxiliary result. 

Lemma 3. Let IX and IX = 1X+A S be quasi- definite vector linear Junctionals and let 7 
and 3 be the generalized Markov functions associated to IX and IX, respectively. The 
sequences of matrix polynomials bi-orthogonal to J and $ are denoted by {KnjmeN; 
{G m } m£ N, {V m }men, and {G m } mS N, respectively. Let also, a™, a™ and /3™, f3™ 
be the leading coefficients of h ofV m , V m and G m , G m , respectively. Then, 

((/cr x c) (acra^Hw 

-K l (0)(/ A r xA r + <54To(z))A T X m+ i(c / ,Q))" 1 ^(^o(z))A T G m (0). (30) 
Proof. Taking into account (|26|) . we have 

V m {h(x)) = D m V m {h(x)) ~ V m {0)5 z (?o{z))A T K m+1 (x, Cl ), 
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where D m — (/3™) /3™ (see theorem 8 of [7]), and thus 
SE(<C) _1 = 2^/ V m (ft(a:))^( a :))G m (/i( a; ))d(ft(i)) 

= f (D m (V m (h{x)) - V m (0) (Inxn +6 z (Mz))k T K m+1 (c l ,c l )y 1 

x6 z (y Q (z))A T K m+1 (x, Cl )))J(h(x))G m (h(x))d(h(x)) 
This means that 

- / M - V m (0) (l NxN + S z (y (z))A T K m+1 (c h q))" 1 5 z (a\,(z))A T G m (0) 
and the statement follows. □ 

Theorem 7. Lei U and U = It + AS be quasi- definite vector linear Junctionals 
and let J and J be the generalized Markov functions associated with 11 and XL, 
respectively. The sequences of matrix polynomials bi- orthogonal to If and 2f are 
denoted by {V m } me w, {G m }meN, {V m }men, and {G„},„ e m, respectively. Let also 
a™, 5™ and /3™, j3™ be the leading coefficients ofV mi V m andG m , G m , respectively. 
Assume that 

lim A m — A, lim B m — B, and lim C m = C, 

where {A m } m ^, {S m } m gN, and {C m } mS N are the sequences of numerical matrices 
involved in the recurrence relation (JTJ) . Then, 

Urn (/Cr^SCKD" 1 =In*n + ?a,b,c(0)(3 < mc (0)- 1 ^b, C (0). 

m—too 

To prove this theorem we need the following lemma. 

Lemma 4. Let {V m } me fj and {G m }meN be a sequences of matrix bi-orthogonal 
polynomials with respect to the generalized Markov function J . Let A m , B mi and 
C m be the matrix coefficients that appear in the recurrence relation §\§ such that 

lim A m = A, lim B m = B, and lim C rn = G, 

tn—too m— too rn— too 

where A and C are non-singular matrices, then 

lim G" 1 ^) = OjvxiV and lim V m 1 (z) = NxN 

m—too m—too 



locally uniformly in compact subsets ofC\T, where T — Hn>oMn, Af/v = U m >ArZ„ 
and Z m is the set of the zeros of V m . 

Proof. From Liouville-Ostrogradski type formula we have that 

— v m n m-l Ur m-2 Ur m-l- 

From the generalized Markov theorem for V m and its analogue for G m , i.e., 
lim V-Hz)*™ ( Z ) = and r lm G« 1 (z)G m 1 (z) = J(z) 

m— too m— too 

we have that 

lim V~ 1 (z)A m 1 _ 1 G m 1 _ 1 (z) =0 NxN . 
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But 

V- 1 (z)A m 1 _ 1 G m 1 _ 1 (z) = V-^(z)V m . 1 (z)V- 1 (z)A- 1 _ 1 G- 1 _ 1 (z) 
According to Theorem [3] 

lim V m - 1 (z)V- 1 (z)A- 1 _ 1 =? AtB<c (z) 

m— >-oo 

and since is bounded for the spectral norm ||.||2 (see Lemma [T]), we have 

that linim^oo G~ x (0) = NxN . To prove lim m ^.oo V^ x (0) = Onxn we follow an 
analogue way. □ 

Proof of Theorem Writing $(q) = (/3™)^ 1 /3™5™(a™)- 1 , by Lemma|3]we 
have 

$m(c;) 

= Inxn - [G- 1 (0)((M,)(3'o(2!)))- 1 y- 1 (0) + G-^O^+^cudW- 1 ^)}- 1 
To analyze lim 3> m (cf) we start by proving that 

m— ¥oo 

lim G-^Oj^iCd.QjJV-^O) = -^B,c(°) J A,B,c(°)^B,c( )- ( 31 ) 

If we put G^(0)K m +x(ci) q))V^ 1 (0) = 7 m (cj), taking into account (J5J and that 
VMz)V-\z) = -VMiV-^z))', then we have 

lm ( Cl ) = A m ^ +1 (0)F m 1 (0) + G m 1 (0)G m+ i(0)C m+1 K 1 (0)(y m 1 (0))' 

= An^x^V-^O) + G m 1 (0)G„ l (0)A„ l K„ + i(0)(F m 1 (0))' 

= -A m [(T/ m (0)^ 1 (0))- 1 (K i (0)F m j 1 (0))- 1 )'(K m (0)^ 1 (0))- 1 ] 

Using dHJ and ([Til), we have (|3 1 [) . 
Finally, we just have proved that 

lim G m 1 (0)((A^)(T (z)))- 1 U m 1 (0) = NxN . 



From Lemma S] we have that lim m ,_ ) . 00 G m (0) = 0at x a, limm^oo K rl 1 (0) = 0at x a, 
and 

||G- 1 (0)((A^)(y o ( 2 )))-^-i(Q)|| 2 

< ]|^ 1 (0)|| 2 ||((A.^)(y (^)))- :l ]| 2 ||V^ 1 (0)]] 2 - 

Thus, the result follows. 

Remark . Notice that in the generalized matrix Nevai class 

lim (/3™)" 1 ^5™(a™)- 



jilt \- J. -jiit in / in y I 

m— >oo 

is well determined and 

lim (/?™)~ 1 ^5™«)- 1 = I™ + S'A,B,c(0)(^ )fl ,c(0))- 1 5'A,fl,c(0). 

// 

5 = /jvxiv + ^,B,c(0)(5 y A) fl, c (0))- 1 5'A,fl,c(0) (32) 
and taking into account that 

lim /?™ = lim (A m ) _1 = lim (C m C m ^ ■ ■ ■ GiA ) _1 = (G • • • GA )-\ 
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lim <= Hm(Q m )- 1 = lim (6 A_i A ■ ■ ■ A m )- X = {Q A ■ ■ • A)" 1 , 
we /iaue that 



lim = lim /3™S lim < = (G- • • GAo^S^^- • • A)" 1 . 

m— yoo m— >oo m— >oo 

Now, we are ready to prove the following asymptotic result. So, we set that 

iim wzr 1 ^ = *■ 

m— >oo 

Theorem 8. Let {V m } me ji be a sequence of matrix orthogonal polynomials with res- 
pect to the generalized Markov function SF and {j4 m } m gN; {5 m } m£ N, and {G m } mS N 
ifte sequences of numerical matrices involved in the recurrence relation ^ such that 

lim A m = A, lim B m = B, and lim C rn = C. 

m— >oo m— >oo m— >oo 

Then, there exists a sequence of matrix polynomials {V m } m ^ left- orthogonal with 
respect to £F defined by (|27l) smc/i i/iai, /or C\{ru {ci, C2, • • ■ , cm}}, we /lave 



lim ^(zJV-^z) 



* _1 [^iVxJV - ^jvxN - 3(3=^^^(0) - (^.B.c^))" 1 )], (33) 



where S zs given by (1321) . 
Proo/. If 

Inxn +S z (7 (z))A T K m+1 (ci,ci), meN, 
is non-singular for every q, a zero of h, then 

F m (/i(*)) = D m V m (h(x)) - V m {Q)5 z {Vv{z))k T K m+1 {x,a), 

where 

tL(O) = D m V m (0)(I NxN +5 z (%(z))A T K m+1 (c h c l ))- 1 

and I? m = Then, multiplying the last relation on the right by V^ 1 we 

obtain 

^(^(a:))^- 1 ^)) = D m {I NxN - V m (Q) (I NxN 

+5 2 (y (z))A T ^ m+1 ( Q ,c ; ))^ 1 ^(a'o(2))A T X m+ i(x,Q))^ m 1 (^(^)) 
But, from the Christoffel-Darboux formula (J7J 

K m+1 (x,ci) = [G m (0)A m y m+1 (/i(a;)) - G m+1 (0)G m+1 V m (Mz))] 
and then, VmV^ 1 can be written as 

V m {h{x))V-\h(x)) = D m (I NxN - 

x (/ A r xj v+^(ToW)A T ^ m+ i(Q,Q)) _1 «5 2 (J , o(2))A T G m (0) 

x (4 m V m+l (Ma:))K n - 1 (/ l (a;)) - G- 1 (0)G m+1 (0)G m+1 ). 
Using (f3U|) in the last relation, we get 

v m {z)v-\z) = [m^fflrHiNxN - \in xN 

+ [(/?™)- 1 ^]K«)- 1 ](A m T/ m+1 (z)t/" 1 (z) - G- 1 (0)G m+1 (0)G m+1 ). (34) 
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If we denote 

E m (z) = G^fOjG^^O)^! - A m 7 m+1 (z))y m 1 (z) 
= C m l +1 G m \MG m {®) ~ V m {z)V m l +l (z)A m l 

then ([33]) becomes 
V m {z)V-\z) 

= IWZr^ZrHlNxN - \lN*N - [{P™y l ~P™][a2{a2r l ]E m {z)). 

From Theorem [3] and its analog for the sequence of matrix polynomials {G m } mg N, 
we have 

lim E m (z) = (Ja.b^O))- 1 - {^A.B.ciz))- 1 . 

m— too 

Since, there exists a sequence of matrix orthogonal polynomials {V m } m £N such that 
lim (/C) _1 /C = we obtain ([33]). □ 

m— >oo 
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